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Reviewing the ideas developed in [lj , the ground state life time of a finite size atomic Bose Einstein 
condensate is studied for coherent, squeezed coherent and thermal coherent ground states. Ground 
state evolution of coherent and squeezed coherent states in a double well potential is studied. Effects 
of thin spectrum on Bose-Einstein condensates is discussed and quasiparticle excitation lifetimes 
are calculated. It is shown that the effect of the states we use on the free energy vanishes in the 
thermodynamic limit. Possible extension to a double well potential and effect of a second broken 
symmetry is also discussed. 
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Shortly after Bose Einstein condensation 0] was ob- 
tained in trapped Alkali atoms [|| , many theoretical and 
experimental studies focused on the quantum coherence 
properties of such systems. It was shown that apart from 
the usual decoherence, which stems from the imperfect 
isolation from the environment, the system also suffers 
phase diffusion [1, [j| , which is due to the atomic num- 
ber fluctuations in the condensate 0]. There is a third 
source of decoherence which limits the life time of excita- 
tions in BECs. This mechanism is based on the existence 
of a group of thin spectrum states Q- The relation of 
thin spectra with decoherence of excitations is recently 
proposed and applied to Lieb-Mattis model and super- 
conductors by the same authors [j| Q . Effect of a thin 
spectrum on quasiparticle excitations on BECs is dis- 
cussed in In this paper we review these ideas with 
some extensions and also discuss the decoherence that a 
double well BEC experiences. 



This paper is organized as follows: Section II begins 
with a review of a toy model for zero mode dynamics, 
studying coherent, squeezed coherent and thermal coher- 
ent ground state lifetimes. A toy model for a double well 
condensate is also introduced and coherent and squeezed 
coherent states in such a system is studied. In section III 
we apply the thin spectrum formalism to atomic BECs 
and discuss its similarities and differences with other cal- 
culations. After making some comments about the exis- 
tence of a second thin spectrum in a double well system 
and outlining the calculation of lifetime in a double well 
system, we conclude in section IV. Acknowledgements 
are in section V. 



In order to understand the basic idea underlying the 
phase diffusion at zero t emp erature, it is useful to in- 
troduce a toy model [l(| The total Hamiltonian of 
a homogeneous Bose Einstein condensate in the weakly 
interacting limit is 
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where a& is the annihilation operator of the k mode and 
u is the parameter determining the strength of interac- 
tions between bosons, that is u — 47rfi ,° s , where V is the 
quantization volume. In order to fix the average number 
of atoms, a chemical potential fj, is also included in single 
particle energy: E k = k^- - /i. 

We are interested in the zero-mode dynamics of the 
system, so we discard terms including &£^ n and redefine 
(X in order to get the basic U(l) gauge symmetric Hamil- 
tonian: 



TL = —a'a^aa — fia^a, 



(2) 



The grounstate of such a Hamiltonian is clearly a Fock 
state with a number N determined by /i. However, Fock 
states have no definite phase and since Bose Einstein 
condensation entails a phase-symmetry broken state, the 
groundstate we seek cannot simply be a Fock state. The 
simplest idea is to consider a coherent state \a) with 
a = y/N, as coherent states are the simplest states car- 
rying a (almost) definite phase. A second step can be 
considering squeezed coherent states [Hj], which again 
carry some phase information, but the uncertainty in 
their phase can be larger or less than a corresponding 
coherent state. And finally, it is important to find a state 
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which both carries phase information and is temperature 
dependent. For this purpose, we are going to introduce 
the thermal coherent states after studying the zero-mode 
evolutions of coherent and squeezed coherent states. 



B. Lifetime for coherent groundstate 

The coherent (or quasi-classical [l3j]) state \a) is de- 
fined as the right eigenstate of the annihilation operator 
and has the Fock state expansion 



harmonic trap frequency uj iT . We find 



t r . = 



N 1 

47T TVcff Wtr ' 
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where N cS = pa\ 



Assuming a typical situation of 



current experiments with N ~ 10 , a s — 10 nm, aho = 
1 pm, and p = 10 21 m~ 3 , we get t c ~ 10/w tr - For a 
harmonic trap with u; tI = 100 Hz, this amounts to t c ~ 
10 _1 seconds, clearly within the regime to be confirmed 
and studied experimentally [||. 



\a) 



(3) 



This simple expansion makes it possible to calculate the 
time dependence of the expectation value of the annihila- 
tion operator a, which we consider as the order parame- 
ter, possible. We define the energy of the n-th Fock state 
|n) as E n = |(n 2 — n) — pn through H\n) = E n \n). This 
gives the simple time dependent expression 
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whose short time behavior is found to be 
(a|o|a) = VNet^e-^e- 1 ^ 1 



(4) 
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Therefore, the order parameter decays exponentially 
0, At longer time scale (a|a|a) revives due to dis- 
crete and thus periodic nature of the exact time evolu- 
tion Q. However, the ratio of the revival time t r to the 
collapse time scales as t r /t c = yN, and hence in the 
thermodynamic limit the collapse is irreversible. Denot- 
ing the density of condensed atom numbers in the quan- 
tization volume as p — N/V, collapse time t c — h/\/~Nu 
can be written as t c — h^/N /puo to see its behavior in 
the thermodynamic limit more directly. As p is fixed in 
the thermodynamic limit where N and V increases in- 
definitely, we see that t c increases indefinitely. Revival 
time t r = h/u = HN/uqp increases with N linearly. In 
practice, the available condensates contain small num- 
ber of atoms and furthermore, they are in traps that 
makes them inhomogeneous systems. Our homogeneous 
system Hamiltonian can qualitatively describe their col- 
lapse time behavior by letting V denote the condensate 
mode volume, though t c ,t r would have different expres- 
sions for the case of a trapped condensate. In particu- 
lar, collapse time of a harmonically trapped three dimen- 
sional isotropic condensate in a coherent state behaves 
like t c ~ TV 1 / 10 in the Thomas-Fermi limit 11]. Homoge- 
nous condensate collapse time is growing much faster, 
as t c ~ N 1 / 2 . Despite these quantitative differences, we 
can still express t c of homogeneous BEC in terms of pa- 
rameters of a trapped BEC. For that aim we shall only 
eliminate m via the characteristic length scale for a har- 
monic trap potential as aho = \J h/mu^ in terms of the 



C. Lifetime for squeezed coherent groundstate 



The squeezed coherent state 
| a, 7) = D{a)S{ r y)\vac) where 



sin) 



12,] |a,7) is defined as 



(7) 



is the unitary squeezing operator and D(a) = exp(cra — 
a*aJ) is the displacement operator. This again is a min- 
imum uncertainty state, but the quantum fluctuations of 
two quadratures are not equal to each other. Arguments 
of 7 and a determine which quadrature is squeezed at 
the expense of increased uncertainty of the other one. In 
particular, if both parameters are real and positive, then 
the state is number squeezed, that is the uncertainty of 
the number operator is reduced whereas the conjugate 
variable, phase, has a higher uncertainty. Such a state 
resembles a Fock state more than a coherent state and 
therefore is expected to have a longer life time, since the 
phase collapse speed is generally proportional to AN, 
which is smaller in this have recently observed 

experimentally [3, EH- This situation is analogous to 
the dispersion of a wavepacket consisting of different fre- 
quency components. 

In terms of a new parameter 



C = 7 



tanh(| 7 |) 
l7l 



(8) 



the Fock state expansion of the squeezed coherent state 
is 



I a, 7) = X^"^'^! 7 
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Here, H n is the n-th Hermite Polynomial. The order 
parameter becomes 



(a, j\a(t)\a, 7) = ^ y/n + lA* n A n+1 e^ E "- 



E n+1 )t 



,(10) 



n=0 



which is not possible to evaluate analytically. We there- 
fore attack the problem using simple numerical methods, 
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and plot the time dependence of the order parameter for 
various values of (, both real and imaginary. It is clearly 
seen in figure[T]that squeezing in the number direction in- 
creases the life time, whereas in figure [5] we see squeezing 
in phase direction leads to a faster decay. 
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FIG. 2: Decay of the order parameter for the coherent state 
and squeezed states of £ = 0.5, £ = 0.5i and £ = —0.5 as 
a function of totr- The same sample parameters with the 
previous plot are used. The solid line is the coherent state, the 
dashed line is the squeezed state with £ = 0.5, and the dotted 
ones are the squeezed states with £ = 0.5i and £ = —0.5. 



FIG. 1: The comparison of the short time decay character for 
a coherent state condensate with that for squeezed states with 
£ = 0.5 and £ = 0.9. Sample parameters used are a a = lOran, 
aho = 1 A 4111 ! n ~ 10 21 rn~ 3 , and a = 10 corresponds to 
N = 100. In this case, the dimensionless time is in units of 
h/u becomes h/u — w^ 1 . The fastest decay (solid line) is for 
the coherent state, while the dashed (dotted) line refers to 
that of a squeezed state with £ = 0.5 (£ = 0.9). 

Considering the time evolution of the Q-functions 
might provide some extra insight into the phase diffu- 
sion. In a contour plot of a Q function, radial distri- 
bution corresponds to the number distribution whereas 
angular one gives the phase information. Beginning with 
a state with some phase information, we expect it to get 
a rotationally symmetric form as time passes and phase 
diffusion occurs. This is seen in the figures [3] and U which 
correspond to squeezed coherent states with a = 10 and 
C = =f0-5 respectively. 
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D. Lifetime for thermal coherent groundstate 

In order to study groundstate lifetime in finite temper- 
ature, we need a state with both a thermal characteristic 
and phase information. The thermal state with the den- 
sity matrix 
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FIG. 3: Time evolution of the Q-function for squeezed- 
coherent state with a = 10 and £ = 0.5 for increasing values 
of ttotr- Figures (a), (b), (c) and (d) correspond to tu)t T = 0, 
tu tr = 0.02, tu tr = 0.10 and tuj tr = 0.40. 



has a uniform phase, that is (a) = 0. We introduce the 
thermal coherent state as p — D(a)p t hD^ (a), where p t h 
is defined using a Hamiltonian without free energy, that is 
E n \n) = (u/2)a ji a jf aa\n). Using the Fock state expansion 
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FIG. 4: Time evolution of the Q-furiction for squeezed- 
coherent state with a — 10 and £ = —0.5 for increasing val- 
ues of tuitr- Figures (a), (b), (c) and (d) again corresponds 
to same time steps; tujtr = 0, tuitr = 0.02, tui tr = 0.10 and 
tuitr = 0.40. The speed of phase diffusion is higher since the 
state is a phase-squeezed state. 
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where L l k are the generalized Laguerre polynomials, 
becomes 
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-^C m (n,a)C* m ,(n,a)\m)(m'\. (13) 



The order parameter, which is found by taking the trace 
of ap equals 




FIG. 5: The short time decays for thermal coherent states. 
The lines correspond respectively to T = 1000 nK, 100 nK, 
10 nK, 1 nK, and 0.001 nK from left to right. The humps 
are due to the ground degeneracy Eo = Ei. Even as the 
temperature approaches zero, the state does not approach 
the ordinary coherent state D(a)\0). Instead, it approaches a 
superposition state D(a)(\0) + |1))/V2- It can be seen that 
the envelope of the function for small T decays at the same 
time scale as a coherent state. 



E. Toy Model for Double Well Potential 

We now consider the case when two condensates in 
identical potential wells are brought into contact via a 
Josephson-like junction. The toy model Hamiltonian is 
of the form M, EDl 



H = -(aVaa- r W&&)-,u(a t a+& t &)-A(a t &+6 t a) (16) 

where a (b) is the annihilation operator for the zero mode 
of the condensate in well A (B). We denote the state 
which has n atoms in well A and m atoms in well B 
by \n,m). Assuming the total number of atoms in both 
wells is fixed and equal to TV, the ground state \gr) of the 
double well condensate can be expanded as 



(a(t)) = V e- 0E -C m (n,a)Cl,{n,a)(k\m)(m'\e% 
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\gr) = ^2c n \n,N - 



(17) 



n=0 



For identical wells, c„ is expected to be peaked around 

/3£ "C m+1 (n, a)C* m {n, a)^^^- 15 ^ . n , ^if/ 2 ' Thc stren S th A of the coupling determines the 

dispersion of the number of atoms in a well. In particular, 
if one makes the ansatz 



There will be destructive interference since E n+ i — E n 
is not constant. The number of contributing factors is 
determined by 0, that is, temperature. Increased tem- 
perature will make more terms contribute and hence lead 
to shorter lifetime. Time evolution of the order parame- 
ter is plotted in figure O for a = 10. 



|c„| 2 oc e (is) 
expanding the Schroedinger equation gives [Tol [Tl1 | 
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FIG. 6: Time evolution of the Q-function of a thermal coher- 
ent state with a = 10 at T = 10 nK. Figures (a), (b), (c) and 
(d) correspond to tu tr — 0, tuj tr = 10~ 4 , t = uj tr = 10~ 3 and 
tutr = 10" 2 . 



| (n(n - 1) + (N - n){N - n - 1)) - fiN. Plotting this 
for coherent, number squeezed and phase squeezed states 
we see that the qualitative results we got from single well 
toy model is still valid, as expected. 
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FIG. 7: Time evolution of the Q-function of a thermal coher- 
ent state with a = 10 at T = 100 nK. Figures (a), (b), (c) 
and (d) correspond to tujtr = 0, tujtr = 10 -4 , tuj tr = 10 -3 and 
tujtr ~ 10 2 ■ Wider radial distribution shows that even when 
the phase information lost, there is a higher uncertainty in 
number due to higher temperature. 



FIG. 8: Comparison of the short time decay for a double well 
condensate which has 200 atoms equally distributed to two 
wells. Solid line corresponds to a coherent state, dashed line 
corresponds to a squeezed state with ( = 0.5 and dotted line 
to a squeezed state with£ = 0.9. 




In order to study the phase collapse, we consider the 
correlation G = (gr\b'a\gr). The time evolution of this 
expectation value after the coupling is turned off can be 
easily shown to be 

oo 

G = J2 e-^ (B "- £ " +l)t <c n+ i A /(n + 1)(JV - n). (20) 

n=0 

We have defined H\n, N — n) — E n \n, N — n), so E n = 



FIG. 9: Comparison of the short time decay for a double well 
condensate which has 200 atoms equally distributed to two 
wells. Solid line corresponds to a coherent state, dashed line 
corresponds to a squeezed state with £ = 0.9 and dotted lines 
to squeezed states with ( = 0.5i and £ = —0.5. 

When the results for single and double well systems are 
compared, it is seen that although decay times are both 
in the same order of magnitude, a double well conden- 
sate with equal number of atoms in each well as a single 
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well condensate suffers faster decay. This agrees with the 
analytical results obtained in [1CJ using the same ansatz 
for single and double well ground states. 



III. THIN SPECTRUM FORMALISM 



He 
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uoPo , t f | 7 \ 



(26) 



Here, uq is the interaction strength not scaled with V, 
that is Uq = Vu, and e k is defined as e k = E k + 2u p . 
The excitation Hamiltonian can be diagonalized to give 



A. Introduction 

By thin spectrum, we refer to a group of states whose 
energy spacings are so low that they are not controllable 
by any experiment and whose effect on the free energy 
becomes zero in the thermodynamic limit. That the ex- 
istence of a thin spectrum leads to decoherence of exci- 
tations at finite temperature is proved in 0, [1[ . In Q it 
is shown that excitations on a Lieb-Mattis system suffer 
decoherence with a rate proportinal to ksT/Nh where 
N is the system size. In [8| the same authors claim that 
this time scale, being independent of the details of the 
system, applies to other physical systems too, and in Q 
they prove that Hubbard model superconductors suffers 
the same decay. In this section, we apply the thin spec- 
trum formalism to atomic BECs and show that they suf- 
fer collapse in the same time scale. 



B. Quasiparticle lifetime in Bose-Einstein 
condensates 



We now go back to the Hamiltonian (fTJ) 



U = E E k a W + 2 E a l+ q a l- q a ka P - (21) 

k k,p,q 



Omitting the 3 rd and 4 th order terms in the non- 
condensed modes (k ^ 0) we get 



H = H Z + H e , 



(22) 



const. 



(27) 



with ui k = [el - UqPo] 1/2 Ell and b k = SakS' 1 . S is the 
multi-mode squeeze operator (l9| . 

For simplicity, we consider a system with quasiparticle 
excitations in only one mode. We denote such a system 
with n atoms in the condensate and m quasiparticles with 
uj by |n, m). Then, 



no|n, m) = n\n,m), 
h k >\n,m) = mSk,k'\n,m), 

H\n,m) = E^\n,m), 



(28) 
(29) 

(30) 



We consider single-particle regime such that quasiparticle 
and particle occupation numbers become the same. Due 
to the number conservation, sum of condensate atoms 
and the quasiparticles should remain the same. To excite 
m quasiparticles, we have to decrease condensate atom 
number by m. The energy of the corresponding state 
becomes 



E {n) = 



2(JV, 



o 



— uqPqti + mu) 



(31) 



Following Q , we assume that in the beginning the sys- 
tem has no quasiparticle excitations at all, and therefore 
has a Boltzmann weighted distribution over the states 
|n,0), i.e., 



p(t = 0)cx JJc 



-P E ° n) \n,0){n,0\ 



(32) 



u 
2< 



H z = ^(nl - h ), 



(23) 



H e = 



k=£0 



(E k + 2uh ) n k + ^ ^a[aL fc a a 



In order to study the excitations and zero mode sepa- 
rately, we need [H z ,H e ] = and for this we neglect the 
quantum nature of clq in H e by replacing tlq/V appearing 
in H e by po = No/V. Here N denotes the number of 
atoms in the zero mode, and so po is the corresponding 
density. After substituting the chemical potential that 
gives the correct number of atoms, uq = uqPq—uoPq/2Nq, 
we get 



H 



U 0P0 -2 - , nj 

— -n - p u n + He 

Z1\q 



(25) 



This state has no phase, and therefore is not the perfect 
starting point for a BEC. However, if we are interested 
only in the collapse of the excitations and if this take 
place on a time scale smaller than the time of phase dif- 
fusion of the zero mode, then we can get an estimate. For 
the time being we will study excitations on this thermal 
state, and in the following subsection we will generalize 
our ideas to thermal coherent states. 

Before proceeding further, we replace the sum in eq. 
(|3"2")l by an integral. Since the value of E n will be ex- 
tremely small for n < 0, it is also legitimate to expand 
this integral to include the negative values of n too. 

Now we bring the system to a superposition of the 
zero-quasiparticle state and the one quasiparticle state, 
that is we bring each \n, 0) to (\n, 0) + \n, l))/^/2. Such 
a state can be interpreted as a particular qubit For 
thin spectrum to affect the system, it is essential that we 
bring the system to a superposition rather than simply 
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exciting a quasiparticle. Now, the off diagonal element of 
the resulting state's density matrix will evolve according 
to 

Pod (t > 0) oc r e-^t'e-^-^dn 

J — CO 

x | e (- fiuopo f2No+itu p /2HNq )n 2 -\-{3p Q u Q n ^ n 

X) 

exp ( r 



,o2 2 2 



2Pu pa/No-2itu p /hN% . 
V \ <~ — . (33) 



^f3u po/2N - itu po/2hN$ 



which gives 



I Pod 



(*)| S 



oc 



tXp \/3 2 N 2 +t 2 /h 2 
+ t 2 /H 2 W 2 



(34) 



after omitting terms with only a phase factor. Although 
the denominator and the numerator have quite different 
forms, we find that both decay in a time proportional to 
t c ~ HNo/ksT. This is the same result that Wezel et. 
al. have found for a crystal [3]. 

For a BEC, we can let N ~ 10 6 - 10 s and T ~ 
l(T 8 -l(r 7 K. This gives t c ~ 10 2 -10 5 seconds, which is 
much larger than even the ground state life times. Unlike 
the room temperature mesoscopic system discussed in 0] , 
BECs are extremely cold systems therefore one single ex- 
citation has such a long life time. However, this does not 
make the calculation unuseful. This life time is given 
only for a single quasi-particle excitation. In general it 
might be useful to have more than one quasiparticles ex- 
cited at a time. If, for example, the condensate will be 
used as the building block of a quantum computer, hav- 
ing an excitation consisting of ~ N quasiparticles will 
make observation of the qubit (superpositions of states 
with m = and m quasiparticles), easier. In order to 
find a decay time for m > 1, the only approximation re- 
quired is 1/N(N — m) ~ 1/N 2 and the timescale will be 
inversely proportional to m. For m ~ Nq, so long as it is 
not the case that 1 — m/N <§; 1, we have t c ~ 10~ 4 — 10~ 3 
seconds. This time scale is much smaller than both the 
observed and expected ground state life times, therefore 
is of interest. 

There are different studies [2(J HH concerning the life 
times of quasiparticle excitations, such as using pertur- 
bation theory, etc. Namely, [2(| has found a linear tem- 
perature dependence for high energy quasiparticles. Our 
calculation will make a quantitative contribution to this 
decay rate. However, low energy excitations are shown to 
have more complex temperature dependencies [2l| . Our 
calculations do not make any predictions for that regime, 
since we have assumed S> uqPq. Also, our theory pre- 
dicts a certain dependence of the life time on the number 
of quasiparticles excited (which is almost linear for small 
m), and this might be used to differentiate it from other 
theories. 

As discussed before, for a group of states to be a thin 
spectrum, their effect on the free energy must vanish. For 



this purpose, we write the partition function as 



Z — Zthin • Z, 



observable 



(35) 



with Z thin = Y in e~P UoPo{n2 '/ 2N °- n \ Again replacing the 
sum by an integral we find that the leading term in the 
free energy per particle \n(Zthin)/N ~ \n(N )/N . This 
means that the mode we consider has no effect on free 
energy and satisfies all the criteria to constitute a thin 
spectrum. 

The next step might be to generalize the calculations 
for coherent of squeezed coherent zero mode occupations. 
However, following this path doesn't give any finite life 
time, since the excitation decay due to thin spectrum 
requires finite temperature, but coherent states have no 
temperature characteristic. But it is natural to do the 
same calculations for thermal coherent zero mode occu- 
pation. Such a calculation is presented in Decay 
rate does not have a linear temperature dependence in 
this case, therefore it might be possible to differentiate 
between thermal and thermal coherent occupations ex- 
perimentally. 



C. Effect of Thin Spectra on a Double Well 
Condensate 

In it is shown that the phase related part of the 
Hamiltonian of a double well condensate can be reduced 
to the form 



U = a+Pi 



a-P 2 



(36) 



Here, a+, a_ and 7_ are parameters depending on sys- 
tem details, P+ and P_ are the momenta corresponding 
to the total and relative phases of the condensates, and 
Q_ is the coordinate corresponding to the relative phase. 
A, again, is a number parameterizing the tunneling be- 
tween different wells. (This Hamiltonian, originally de- 
rived for a condensate consisting of two different types of 
atoms in a single potential well, is applicable to a double 
well system when the parameter corresponding to col- 
lisions between different types of atoms is taken to be 
zero.) When coupling A between the wells is taken to be 
zero, this Hamiltonian reduces to one of two free parti- 
cles: 



U = a+Pj_ 



a- Pi 



(37) 



It is seen that the system has two modes corresponding 
to motions without restoring forces. The reason is that 
now that there are two wells, there is an extra symmetry 
that is spontaneously broken. Therefore, there are two 
different thin spectra. 

Cumulative effect of multiple broken symmetries (and 
hence multiple thin spectra) on excitation lifetime is 
studied in [l[ . If lifetime corresponding to individual thin 
spectra are t\ and £2, then the resultant collapse time is 
the harmonic sum of individual life times: 



to 



(38) 
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Using a's corresponding to the system under consider- 
ation, it is thus possible to find the resultant life time, 
which is supposed to be in the same order of magnitude 
with the smaller life time. 



IV. CONCLUSIONS 

Generalizing the Toy model calculations [1(1 fllj . we 
discussed the phase decoherence of coherent, squeezed 
coherent and thermal coherent states. For visual clarity, 
time dependence of various Q functions is shown. A gen- 
eralization of the toy model to double well systems is also 
discussed and time evolution of the order parameter is 
studied for coherent and squeezed coherent states. This 
step, being important not only for double well BECs, 
might bear important results for any Josephson-coupled 
system. 



The effect of thin spectrum [7|, l8| on quasiparticle ex- 
citations in BECs is briefly reviewed. It is shown that 
the presence of the so called thin spectrum states, which 
have vanishing level spacing, also has no effect on free 
energy per particle in the thermodynamic limit. Qualita- 
tive dependence of life time on the number of excitations 
is given. Finally, as a simple example of a system with 
more than one spontaneously broken symmetry [l| , a cal- 
culation of excitation life time in a double well system is 
outlined. 
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